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Abstract 

In-medium interactions of a particle in a hot plasma are considered in the 
framework of thermal field theory. The formalism to calculate gauge invariant 
rates for photon and dilepton production from the medium is given. In the 
application to a QED plasma, astrophysical consequences are pointed out. 
The photon production rate from strongly interacting quarks in the quark- 
gluon plasma, which might be formed in ultrarelativistic heavy ion collisions, 
is calculated in the previously unaccessible regime of photon energies of the 
order of the plasma temperature. For temperatures below the chiral phase 
transition, an effective field theory incorporating dynamical chiral symme- 
try breaking is employed, and perturbative QCD at higher temperatures. A 
smooth transition between both regions is obtained. The relevance to the soft 
photon problem and to high energy heavy ion experiments is discussed. 
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I. INTRODUCTION 



Considerable effort is invested in present and future experiments of ultrarelativistic heavy 
ion collisions (URHIC) in order to observe an excursion of the bulk of strongly interacting 
matter from the state of hadrons before the collision into the phase of a quark-gluon plasma 
(QGP) [|IJ. In order to see this shortlived state directly, one would like to observe photons 
emitted from the hot plasma, as well as dileptons. Since these probes interact only electro- 
magnetically, their signal is not distorted by later interactions as are other particles which 
are studied for the same purpose. 

The experimental capability of measuring electromagnetic probes was demonstrated in 
the photon channel by Helios @, WA80/98 @ and CERES Q as well as in the dilepton 
channel by Helios § and CERES §. 

The signal originating from the plasma phase is, however, buried under a background of 
photons from different origin such as from the decays of 7r° or rj or from hadronic reactions 
at a temperature comparable to that of the deconfmed phase ]7). After subtraction of these 
sources, a remaining signal seems to persist in part of the experimental analysis. At present, 
it is vividly discussed to what extent one can account for these data within more H or less 
|| conventional physical pictures. 

However, our theoretical knowledge of the spectrum of electromagnetic probes from both 
the plasma as well as the hadronic phase is still uncertain to some extent. A better handle 
on these spectra from theoretical calculations is necessary in order to disentangle the various 
sources and to identify the phases reached during the collision. In particular for the soft 
part of electromagnetic radiation, this problem represents a challenge to theory in itself, 
due to the nonperturbative nature of the photon emission process: Multiple rescattering of 
the emitting particles and the Landau- Pomeranchuk-Migdal (LPM) effect play an important 



role in the medium for photon energies E 1 <T |TQ , |TT|| , as well as for dileptons of an invariant 
mass in this range. 

This problem motivated the present work, in which we will investigate the production 
of photons and dileptons from a strongly interacting plasma at finite temperature. After 
a short sketch of the insufficiencies of existing calculations, we show how one can reach an 
improvement by taking thermal scattering and thereby spectral broadening of the emitting 
particles in the heat bath into account. The problem is addressed in the framework of 
thermal field theory, results are given for a QED plasma as well as for a QGP within a 
model incorporating dynamical chiral symmetry breaking. A part of the results has been 
presented already in a short paper |12| . 



In an a s expansion, the lowest order of photon production proceeds via annihilation 
(qq — > g-y) and Compton (qg — > qj) processes. In next to leading order (NLO), numerous 
corrections to these processes arise, a complete calculation of the order 0(aa1) has been 
achieved in [□[]. With an initial quark and gluon distribution specified by distribution 
functions fi(Ei) and / 2 (£ , 2 ), and the final state quark or gluons distribution f 3 (E 3 ), the 
production rate reads 

dN° 
d 6 p 
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where the last factor takes into account Pauli blocking or Bose enhancement of the quark 
or gluon in the final state, and M stands for the elementary cross section considered. 

The production of hard (high p±) photons in reactions of colliding hadrons has been 
calculated using 'cold' parton distributions delivered by structure functions and using M in 
NLO. For sufficiently high p±, very good agreement with the corresponding data is reached 
fl4| , only towards low p± some discrepancy has been reported ||15|| . This may hint at the 
insufficiency of using even NLO calculations in the soft regime, but may equally well be due 
to our still insufficient knowledge of the parton distribution functions in the relevant x and 
Q 2 range, see |T(| for an analysis. Even with this minor uncertainty, one has reached a very 
good quantitative understanding of photon production. 

Now let us look at the same processes in a plasma, where the partons have reached a 
thermal distribution. We will consider situations in which the spatial extension of the plasma 
is lower than the mean free path of the photons emitted, i.e. we consider the emission of 
'white' radiation in contrast to thermal black body radiation. Due to the small size of nuclei 
compared to the mean free path of an electromagnetically interacting particle, this is always 
the case for heavy ion collisions. Using M in lowest order, and taking thermal quark (q), q 
and gluon distributions of temperature T results in a production rate R (per unit volume 



element) as [17 



E dN ° 



5 aa, 



d 3 p 9 2tt 2 



T 2 e -E/T 



log 



ET 



77T 



(2) 



with some constant c°. This rate diverges when m — > 0, which is the crucial limit of 
chiral symmetry restoration for strongly interacting quarks approaching the phase transition 
temperature. This unphysical divergence will eventually be shielded by medium effects on 
the emission process. 

A step towards the calculation of such medium effects has been the application of the 
Braaten-Pisarski method of hard thermal loops |Tj| to this problem [0,|T^|. The resulting 
photon production rate is 



5 aa s ~ p/r , C\E 



(3) 



with a constant c\ ~ 3 and the strong coupling constant g. For g 2 ~ ci, the term log(E/T) 
reminds us of the validity of this approach only in the region of E 1 ^> T. 

A more detailed investigation of this infrared problem within the hard thermal loop 
(HTL) method of Braaten and Pisarski has been presented in In these works, the 



production of soft photons was studied thoroughly. The origin of the infrared problem could 
be traced back to divergences which occur when the real photon is emitted collinear to 
a thermal gluon. For this reason, it was concluded in |^U|,^T| that no finite value for the 
production rate of soft photons can be obtained within the HTL method. 

This is the motivation for the present work. Physically, the HTL method takes into 
account the thermal masses particles acquire in the medium. In addition to that, we now 
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also consider the thermal scattering of the partons, which results in an energy uncertainty 
as the quark propagates. As we will show below, this is the dominant physical process 
to be considered for quarks emitting soft photons with energies E 1 < T. The thermal 
scattering leads to a finite lifetime (or nonzero spectral width) of every excitation in the 



medium, described in analogy to the decay width of an excited state ||22|| . One effect of 
such a spectral width is that it naturally removes the infrared divergences mentioned before, 
therefore enables us to calculate production rates for soft photons. Secondly, the energy 
"uncertainty" , which is actually a kind of Brownian motion, is directly related to the emission 
rate of thermal photons. 

This paper is organized as follows. In the next section, we present the general formalism 
of thermally scattered particles in a heat bath. This includes the calculation of thermal 
widths as well as the photon production rate in a gauge invariant manner. Section II. C 
gives the comprehensive example of a fermion (quark) in a QED plasma. Although not 
being realistic for the QCD case, it allows for simpler and often analytic solutions and thus 
for a clear illustration of the relevant physics. 

We then turn to the case of the QCD plasma, section III. A crucial aspect of QCD 
at low temperatures is the breaking of chiral symmetry. Hence up to the chiral phase 
transition temperature we describe the plasma by the Nambu-Jona-Lasinio model, which 
incorporates this feature dynamically (Sect. III. A). In the subsequent part of this work, 
III B, we turn to high-temperature perturbative QCD, with temperature dependent strong 
coupling constant a s . Section III.C gives the results for the photon production rates over 
the entire range of temperatures and of photon energies, and we discuss the relevance to a 
variety of experimental situations in section III.D. 



II. PHOTON RADIATION IN THERMAL FIELD THEORY 

In this section, we first briefly recall the formalism of thermal field theory using spectral 
functions, outline how the self energy of a thermal particle is obtained in general, how it 
is related to the thermal width and how gauge invariant rates for photon production are 
obtained therefrom. We finally illustrate the achievements with the example of a fermion in 
a QED plasma. 



A. Spectral functions and self energies 

For any physical system one would like to have a causal description: Physical particles 
e.g. may exert a measurable influence only after their emission. In the framework of quantum 
field theory this means that one would like to use causal Green functions or propagators 
in the theoretical description. The requirement of causality however touches two aspects 
of field theory. It relates the boundary condition in time that a propagator fulfills to the 
average occupation number of the state that is propagated. 

For a vacuum state, this leads to the well-known Feynman boundary conditions, which 
in terms of the free propagator in momentum state translate into the simple +ie- description 
in the denominator. 
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At nonzero temperature, the average occupation number of a state is given by a ther- 
mal equilibrium distribution function (Bose- Einstein or Fermi-Dirac). Hence, the temporal 
boundary conditions for the propagation of particles at finite temperature are more compli- 
cated than in a vacuum state, they are called the Kubo-Martin-Schwinger (KMS) condition 

This KMS condition leads to a causal propagator with a complicated analytical structure. 
It is therefore safer for thermal systems to deviate from the description in terms of causal 
propagators. Rather one uses only retarded and advanced propagators, whose temporal 
boundary conditions do not depend on the occupation number of states. It is well known, 
how to express a finite temperature perturbation theory in terms of retarded and advanced 
propagators (see refs. |2^ , |25[| for an extensive discussion). 



Mathematically, the retarded and advanced propagators are analytical functions of their 
energy parameter in the upper or lower complex half-plane. Analytical functions however 
obey the Kramers-Kronig relation, and this implies that the retarded propagator of an 
interacting field theory is known completely if only its imaginary part (or spectral function) 
A is known along the real axis. Hence, for the retarded quark propagator in our system we 
write, for arbitrary complex energy E 

oo 

S RA (E, p) = J alE> A q (E', p) E _l l±ie ■ (4) 

— oo 

For free particles the spectral function is proportional to a 5-function, 

4 ree = (£ 7 ° - PI + m q ) sign(E) 5{E 2 - p 2 - m 2 q ) , (5) 

which kinematically limits the asymptotic states to be on-shell. 

However, as we have argued above, such asymptotically stable states are not present in a 
finite temperature system: Each particle is subject to collisions which will add a statistical 
(thermal) uncertainty to its energy as function of time (thermal scattering, or Brownian mo- 
tion). This indicates that the limit of a 5-like spectral function cannot be used in interacting 
thermal systems - which has been proven rigorously in the Narnhofer-Thirring theorem p2| . 



One may attribute this to a fundamental property of temperature: A thermal particle 
distribution function has a special rest frame, hence corresponds to a breaking of the Lorentz 
invariance. It is well-known that a state of broken symmetry requires to chose adequate 
basis functions for a quantization. In case of the finite temperature breakdown of Lorentz 



invariance, the basis functions turn out to be quantum fields without a mass shell p6| , |27 
In other words, the field theoretically correct way to treat a finite temperature system is in 
terms of continuous spectral functions. 

How to put these two aspects together, i.e., the transformation to retarded/advanced 
propagation as well as the perturbative expansion in terms of generalized free fields with 
continuous mass spectrum, is discussed in ref. [^]. For the purpose of the present paper, it 



is sufficient to choose a parametrization for such a spectral function. 
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B. Parametrization of the quark spectral function 



For this parametrization we take as a guideline the idea to be not too far from the 
quasi-particle picture, i.e., we make an ansatz for the inverse retarded quark propagator 

P^Y -m°- Z R {p) « (po ± i7g)7 - p~y - m q (6) 

with a given self energy function Ti R (p) in the vicinity of po = ±^/p 2 + m 2 and \p\ <C m 



This implies, that we assume the whole model to be dominated by its infrared sector, see 
the remark at the end of this subsection. This ansatz translates into a spectral function as 

= 7,7o(^ + ^(p) 2 )-2^ P + 2^ 
9 ?r (£ 2 -^(p) 2 ) 2 + 4£ 2 7 2 

Here, 7 M = (70, 7) is the four- vector of Dirac matrices and £l q {p) 2 = p 2 + m 2 + 7 2 . m q is the 
dynamical mass of the quark, and its spectral width parameter we label j q . Note however, 
that the half-maximum width of the spectral function peak is 2 , ~f q . 

One may regard this spectral function as the generalization of the standard energy- 
momentum relation of eq. (||) to a broader distribution for thermally scattered particles, in 
this particular case represented by a double Lorentzian. Note also, that this parametrization 
differs from a quasi-particle approximation only by one parameter 7^, and in the limit 7 g — > 
one recovers the free spectral function @. 

Moreover, it may be shown explicitly, that this spectral function has a four- dimensional 
Fourier transform that vanishes for spacelike coordinate arguments [ 28fl . Since the Fourier 



transform of the spectral function is nothing but the expectation value of the anti- 
commutator function of two quark fields, this is an important aspect: It guarantees, that 
there cannot be any propagation of interactions faster than light. Let us note, that with a 
general momentum dependence of the spectral width parameter, this requirement may be 
violated. 

For the self energy function we use expressions obtained in a skeleton expansion of the 
full Green function, i.e., we employ Feynman diagrams for this self energy which are again 
functionals of the spectral function we wish to determine. 

In such an expansion, the one-loop (Fock) diagram, depicted in fig. [j], is the lowest 
order term with a non- vanishing imaginary part. In the following, we restrict ourselves 
to this lowest order. We consider a model where quarks are coupled to different types of 
bosons, to be specified later. The calculation of the Fock self energy with full propagators 
is straightforward |25| and gives for the imaginary part 

lmZ R (po,p)= (8) 

00 

-7T Jj^y, JdET fl A q (E,k)T u A B u (E- Po ,k-p) (n q (E) + n B (E - p )) . 
—00 

Here, Ab is the boson spectral function, V ^ and T u are the interaction matrices at the 
vertices, and ns {n q ) is the standard thermal equilibrium Bose (Fermi) distribution functions 
at temperature T, 



6 



P = (Po,P) 



£1 




FIG. 1. Self energy diagrams for the photon production rate. 
Left: Fock diagram for the quark self energy contribution, 

right: photon polarization diagram including full fermion propagators on both lines. 



*B, q (E) = -pz— . (9) 
e -\- i 

The real part of this self energy function is determined by a dispersion integral, similar to (|J) 
for the propagator. Note, that the divergence of this integral either requires renormalization 
or a regularization procedure. 

Having specified the self-consistency criterion for the quark propagator, we may now 
ask for its validity. In particular, one may suspect that representing the complicated quark 
spectral function over the whole range of energies and momenta by only two parameters 
is an oversimplification. However, we find on the contrary that for the self-consistent Fock 
approximation with massless vector bosons the fermion spectral width is dominated by 
a constant term [ 25| . Hence, at low temperatures our ansatz for the spectral function is 
consistent for quark momenta \p\ < m q . 

For higher temperatures the quark mass is small, whereas the quark momenta are typi- 
cally of the order of the temperature. However, as has been shown in ref. PD| , the damping 



rate is dominated by the minimal distance in the complex energy plane between the origin 
and the spectral function pole (note, that physical propagators do not have poles in the 
complex plane). This minimal distance is again given by the width, see eq. (|29|) . 

For the loop integrals in self energy functions the limitation to small quark momenta is 
in principle violated. However, equilibrium distribution functions effectively provide a cutoff 
at momenta |p| ~ T. We therefore find, and have confirmed this by extensive numerical 
computations, that the ansatz of a momentum independent spectral width parameter is 
very well justified for temperatures T ~ s/n^g + 1\ ~~ a relation, which is satisfied in our 
approach. Only in the limit of asymptotic freedom, where the coupling parameters indeed 
become small, this approximation will possibly fail. 

C. Gauge invariant photon production rates 

The width calculated from the quark self energy diagram now enters the photon polar- 
ization n at finite temperature, see fig. [L| The imaginary part of the retarded one-loop 
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polarization function U R is p5 



Imn^(fc ,fc) = -7reJ / / dE 



— oo 



Tr^A^ + ^P + fch.AO^P)] K(^)-n,(E + fc )) , (10) 

where e q is the electric charge of the quark. The photon production rate for the hot plasma is 
proportional to this imaginary part, summed over the different physical photon polarization 
directions. 

We now address the question of gauge invariance of the rate calculated in this manner. 
The photon production rate is gauge invariant if the current which produces the photons 
is conserved. For the current conserving polarization tensor, which we denote by EE, this 
implies transversality, k^T\.^ v {k) = 0. 

The polarization tensor II as calculated from eq. (|l0|), which is connected to the current- 
current correlator U tJiU (x, y) oc (VvfyV'a: ■ ipyjutpy), violates this requirement, and may not in 
general be used for a calculation of the photon production rate. 

This can be traced back to the fact, that the naive current ipj^ip is not conserved. Of 
course, a theory with a nontrivial spectral function also has a conserved (electromagnetic) 
current - but this differs from the naive expression |31|. 



Let us briefly discuss the nature of this difference, starting from the lagrangian of a 
generalized free field which gives rise to a propagator with certain self energy insertion. A 
detailed discussion is carried out in ref . |TIJ . For a one-component fermion field, this would 



be 

£[V>]=^(x)(i^-mo)V>(z)- Jd 4 y^P(x)J:(x,y)^(y) . (11) 
Performing a local phase transformation of this field then allows to find a conserved current 
f(x) =i){x)^ip{x) -i fd A yd 4 z^{z)A^{z,y-x)^{y). (12) 



The function A M is a vertex correction function. Current conservation now is equivalent to 
the fulfillment of the Ward-Takahashi identity, which for the Fourier transformed quantities 
reads 



{p-q)^(p,q) = Z{p)-Z{q) . (13) 

In this equation, p — q = k is the photon four momentum. Without loss of generality we 
may fix the photon 3-momentum to be the vector (0, 0, k). It is then obvious, that the Ward 
identity involves only the components A and A 3 , it does not restrict the transverse parts 
A 1 and A 2 of the vertex correction function. In short words, the logical steps are: 
A , A 3 nontrivial =>- transversality of the polarization tensor II =>- current conservation =>- 
gauge invariance of the photon production rate. 

Naturally this does not imply, that vertex corrections - if calculated diagrammatically - 
have only A , A 3 components: It merely tells us, what is sufficient to ensure gauge invariance 



S 



of the photon production rate. Specifically for the spectral function we have postulated, only 
A is necessary to acquire a conserved current. 

Correspondingly only the components U 0u = U u0 of the current conserving polarization 
tensor are different from the one-loop expression fllPf) . This tensor is the autocorrelation 
function oc (j^ j u ) of the conserved current. It is crucial to realize that the space-like 
components are not modified, IP- 7 = IP- 7 . 

In the next step we use this fact together with the condition of on-shell photons, ko = \k\. 
Current conservation implies that k^Uoo — | | 2 II33, i.e., Hoo = II33 for on-shell photons. 
Hence, these two components cancel in the sum over polarizations: 

Let us note, that this chain of arguments is rigorous: One may debate, whether vertex 
corrections are necessary out of phenomenological reasons - but they are not necessary in 
order to achieve a gauge invariant result for the on-shell photon production rate. Of course, 
this conjecture has a drawback: For the more general case of a momentum dependent spectral 
function the argument above does not hold any more. Also the calculation of off-shell photon 
production, such as required for dilepton production rates, necessitates the calculation of 
vertex corrections. 

To summarize this discussion: On the level of our approximate spectral function, and 
with the fully causal propagators following from this spectral function, the gauge invariant 
photon emission rate out of the hot plasma is 

r(e„ t) = = 2 n -^fl im (n« + n* ) = JL (n< + n<) , (is) 

where Il < stands for the off-diagonal component in the standard notation of real-time ther- 
mal field theory ||25|| . 

The polarization tensor U(x, y) is the correlator of electromagnetic currents at different 
space-time points. Hence, interference effects between photons emitted from different points 
in space and time, as far as they affect the single photon rate, are automatically taken 
into account. This also includes the Landau- Pomeranchuk-Migdal (LPM) effect: Multiple 
thermal scattering of the slowly moving quarks leads to the interference of sequentially 
emitted soft photons, thereby reducing the soft photon rate. The equivalence of the semi- 
classical LPM description with the field theoretical formulation used here has been proven 



in ref. [11 . We will show that it is exactly this interference effect which gives rise to our 



primary result. 



D. QED plasma and astrophysical consequences 

In this subsection, we illustrate the physics with the example of a fermion (quark, or 
"heavy electron") interacting electromagnetically in a plasma of temperature T. This can 
be handled to a large extent analytically and thus allows for a clearer understanding of the 
mechanism we are treating. 

For the purpose of this example, we keep in the following the mass of the fermion fixed 
at an (arbitrary) value of 300 MeV, and solve the self consistent equation of the in medium 
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FIG. 2. Photon production rate i? 7 from an electromagnetically interacting particle of 300 MeV 
mass in a plasma as a function of the photon energy for different temperatures T. 

quark propagator, eq.(|6]), for the width only, using using real and imaginary part of the Fock 
self energy (§) for massless free photons. The cutoff for the determination of the real part of 
the self energy function is chosen as A + a/A 2 + ra|, with A = 650 MeV (see the discussion 
in section III. A). 

Although this width 7 em is, in principle, a non-analytical function of the temperature 
29fl , the smallness of the electromagnetic coupling constant allows to approximate it very 



well by the lowest order result 

7 ^(T) « ^ . [x _ W 5 x _ 10 a AT 5 

lq v ; 9 1 m q J 9 L 9 vr m\ J 9 V ; 

where ReS y is the Lorentz vector (oc 7 ) component of the fermion self energy function, 
and the factor 5/9 is due to the (u,d)- family averaging of electric charge. The purely 
electro dynamic quark spectral width as function of temperature is plotted in fig. |], together 
with the other contributions to be discussed later. 

The photon production rate we obtain from eq. ( |T5| ) with this width is shown in fig. |2] for 
various typical values of the temperature. For small photon energies, i.e. very soft photons, 
we find a saturation of the rate below values of E 1 = 2 / -f^ m w 2 ■ 5/9 • aT. The factor 2 
arises, because the half-maximum width of the Lorentzian spectral function peak is 27^ in 
our choice of parameters. 

The physical interpretation of this effect is obvious: The emission of low-energy photons 
requires unperturbed propagation of the emitter over the wavelength of the photon. Along 
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its path however the quark is subject to thermal perturbations - and this hinders the photon 
emission for for E 7 < 2 r y q . Our result agrees qualitatively with the conjecture of Weldon [10 



Moreover, we could clarify the dominant suppression scale to be twice the spectral width 
parameter of the emitting particle, or equal to the half-maximum width of the spectral 
function peak. In the spirit of the last remark in section II. B, we may state that this is an 
interference effect between photons from different points in space and time. 

The rate for high energy photons falls off with photon energy E y as e~ E ^l T from the 
Boltzmann factor. For a particle mass m q ^> 7 g this result coincides with previous calcu- 
lations of eq.(|3|) with equivalent parameters. However, in contrast to this calculation, our 
result does not employ singular behavior in the limit of m q — > 0. This illustrates nicely how 
the finite thermal particle width regulates the infrared behavior. 

The photon production rate may be approximated as 

Rl = ^P 3 —, e-*^ z[T] , z[T] oc ( T f 1 f ^ . (17) 
fit E 2 + 4 7 2 \ T for E 1 > 2 7g v ; 

For all temperatures, the limit E y — > oo is determined by the Boltzmann factor e~ El ^ T . 
Note, that this functional dependence is a fitted result after the self-consistent calculation, 
hence 7 g is not an external parameter that my be varied independently of T and the coupling 
strength. Numerical fits will be published separately. 

Another remark may be appropriate as to compare this result with standard (semi- 
classical) treatments of the Landau- Pomeranchuk-Migdal effect: Although the exact value 
of 7 g may be debated, the general form for the radiation rate we obtain conforms very well 



with all the semi-classical approaches [11 



Let us briefly consider the astrophysical consequences of this result. For this we regard an 
era of the universe, where the photon energy density dominates the matter energy density, 
roughly 1 MeV < T < 10 MeV. Using eq. (fL7j) together with eq. (p~6|) , we find that the 
photon number density in the early universe is roughly n 7 Pd 2/tt 2 ■ T 3 , while the photon 
energy density is given by p 7 ~ 6/n 2 ■ T 4 . These have to be compared to the standard 
Bose-Einstein values of n® « 2A/n 2 ■ T 3 , p 7 « 6.49/7T 2 • T 4 . This comparison implies, that 
the effective number of degrees of freedom, g* in the total energy density 

^P T ' < 18 > 

is somewhat reduced by the coherence effects we are considering. This then would lead to 
a faster expansion of the universe during the time where our result applies. However, this 
effect is on the order of a few percent and therefore in the moment beyond the reach of 
experimental observation. 

Furthermore we may assume, that in the present spectrum of the cosmic background 
radiation the photons have retained the spectral cutoff point from the moment when the 
universe became transparent. If we assume this transition to happen at a temperature of ~ 
1 eV, the effect we are proposing predicts the cosmic microwave background to be thermal 
at wavelengths shorter than A cu t = rr/aT, but leads to a reduction of the long- wavelength 
background radiation below its black-body value for larger wavelengths. Due to the expan- 
sion factor of ~ 1000, we estimate the present cutoff wavelength to be approximately A cu ^ ~ 
1 m, which makes it difficult to observe this deviation. 
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FIG. 3. Diagrammatic equation for the quark propagator at low temperatures. 
Double line = effective quark propagator, wavy line = photon. 



III. PHOTON PRODUCTION FROM A STRONGLY INTERACTING PLASMA 

We now come to the main topic of this work, the production of photons from a QGP. For 
this purpose, we distinguish two temperature regimes. In the low T region, below and around 
the phase transition temperature, dynamical chiral symmetry breaking and its restoration at 
T c ~ 200 MeV plays an important role and has to be incorporated in a realistic description 
of the quark dynamics. We do so by using the Nambu-Jona-Lasinio model as an effective 
model up to T ~ T c , and add to this model the self-consistent summation of quark-photon 
Fock diagrams. 

In the high temperature limit, chiral symmetry is restored and the coupling is small 
enough for a perturbative expansion in the strong coupling constant. In this region, we 
therefore use a self-consistent determination of the quark width obtained in perturbative 
QCD0. 

In both regimes, nontrivial spectral functions for the respective interacting bosons are 
used. 



A. Nonperturbative temperature regime 



Here, we consider the Nambu-Jona-Lasinio model |$3| in the SU(2) version on the quark 



level, see |35| for a review and the notations used in the following. 

In this effective field theory, the nonperturbative interaction between quark and antiquark 
fields at low momentum transfer is described by the sum of a scalar and a pseudoscalar local 
interaction, C int = G[(^ (x)^/ (x)) 2 + (^(x)i , y 5 T'^(x)) 2 } where \I/ = (u,d). This is understood 
to be a summation of nonperturbative gluon interactions among the quark fields. It models 
the chiral symmetry properties of QCD in the nonperturbative regime, which is essential to 
address processes on the energy scale of the temperature. 

This effective field theory models the chiral symmetry properties of QCD in the nonper- 
turbative regime by a quartic self-interaction of quarks. At small temperature, the dominant 
contribution to the quark self energy is the tadpole (Hartree) term, which is expressed in 



terms of the spectral function as [[25[] 



E H = -2GN c N f j-^- JdETr[A q (E,p)}n q (E). (19) 
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FIG. 4. Contributions to the width 7 ? of a quark. 
Short dashed line = scalar meson, long dashed line = pion. 



Like any nonrenormalizable model, the NJL requires a momentum cutoff A, which can be 
motivated as a crude incorporation of asymptotic freedom at large Q 2 . For the present 
generalization, this cutoff is shifted to the energy integration, 



A q = y /A? + m^. (20) 

Usually, the temperature dependent quark mass m q (T) is the solution of the gap equation 
m q = m + T, H \m q ) . With appropriate parameters, this describes the scenario of spontaneous 
chiral symmetry breaking, i.e., the transition from a current quark mass mo ~ 5 MeV to the 
constituent quark mass m q w 1/3 x the nucleon mass, and its restoration at a transition 
temperature T c . The only parameters were chosen as tuq = 5 MeV, A = 0.65 GeV and 
G = 5.1 GeV -2 , and result in T c = 193 MeV and a vacuum mass of = 140 MeV for 
the pion and m q = 332 MeV for the constituent quark (135 MeV resp. 331 MeV without 
photons) . 

The Fock self energy is the next-to- leading order contribution in a 1/N C expansion [|36 



We consider quarks with four-momentum (po,p) = (m g ,0), hence we can decompose the 
retarded Fock contribution to the self energy in a (complex) scalar and a vector part as 

s Fock = V jS + 7oE V. 

The photon Fock contribution is added to the Hartree self energy, and instead of the gap 
equation we solve eq. @ for the mass and width of the effective quark field. Split into real 
and imaginary part this reads 

rrtq^iq = — Im£ (m q , 0) (m q — ReS y (m 9 , 0)) 

- ImS 5 ^, 0) (mo + S H + ReZ s (m q , 0)) 

-7? = (m, - ReZ v (m q , 0)) 2 - (m + Y> H + ReE 5 (m„ 0)) 2 

- ((ImS y K, 0)) 2 - (ImE 5 K, 0)) 2 ) . (21) 

Here, ImS v is an even function of po, ImE 5 is an odd function of p^. For 7 g — > 0, the above 
equations reduce to the standard gap equation of the NJL model. 

The NJL model furthermore describes two bound states which are regarded as effective, 
T-dependent pseudoscalar 7r-meson and scalar a-meson. 

The traditional way to determine the masses of these collective modes is to solve the 
equation p5j 
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1 -2GRen(m B ,0) = , 



(22) 



with II the corresponding polarization function determined as a dispersion integral over eq. 
( [10D . As a cutoff for the dispersion integral one uses ±2A g according to eq. (|20|) . 

Within our formalism we use the following ansatz for the retarded boson propagator 
along the real axis: 

- m 2 - Tl R (k) = {k - (u B (k) - i 7B )) fa + (u B (k) + i 7B )) (23) 

with uj 2 b = m 2 B + k 2 . This translates into a spectral function as 

A B {E,k) = - -, (24) 

1 ; vr (e 2 - n B (ky) 2 + ae^i 1 ; 

with Q 2 k = u B (k) 2 + 7 |. 

The parameters m B and 7b are determined from the equations 

1 - 2G Ren(m B , 0) + (G vrcx(m B , 0)) 2 = 

m B Gna(m B ,0) = j B . (25) 

The mesonic Fock contributions to the quark self energy are treated only perturbatively, 
i.e., their imaginary parts are used to modify the quark width according to fig. |4] and their 
influence on the quark mass is neglected B9j. 

To check the consistency of this approximation, we also performed a self consistent sum- 
mation of the meson Fock diagrams, which gives rise to a small correction of the constituent 
quark mass as well as the meson masses in our extended NJL model. However, such a 
summation breaks chiral symmetry explicitly - and one may expect, that the mass shift is 
canceled by other diagrams. For the purpose of the present paper we therefore chose the 
perturbative treatment of the meson Fock self energy contributions as described above. 

Physically, our approach amounts to consider photon emission processes, which are ini- 
tiated by the interaction of the quark with a single hot meson. The resulting quark width 
7g is plotted in fig. For low temperatures, we again find 7 cx T for each channel. Up to a 
temperature of ~ 150 MeV, the quark width 7g is in fact dominated by the purely electro- 
dynamic contribution. This indicates, that photons should be taken into account even for 
strongly interacting systems at such temperatures. 

Due to the quasi-Goldstone mode of the pion, its contribution to the quark width remains 
negligible up to the Mott temperature Tm = 212 MeV, which is defined by m n (TM) = 
2m q {TM) as the point where the pion can dissociate in a qq pair. For T > Tm, the pionic 
contribution to the quark spectral width is actually dominant. Towards higher temperatures, 
the competing effect of an increase of the mass of the 7r (now a resonance) again turns the 
width down. 

One may argue at this point, that in the NJL model quarks are not confined. However, 
the above results may be translated to other models as well: They represent nothing but 
a critical opalescence to photons in the vicinity of the chiral phase transition. Hence we 
expect the drastic increase of the effective 7 around T c to be quite independent of the model 
used. 
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FIG. 5. Feynman diagrams used in the calculation of the effective quark and gluon propagator. 
Each line represents an effective propagator here, dashed are the ghost fields. 

B. Perturbative QCD for high temperature plasma 

In the high temperature limit, a calculation within perturbative QCD becomes sensible. 
Furthermore, non-abelian gauge invariance becomes an imperative of the calculation. For 
each of the degrees of freedom, i.e., quarks, transverse and longitudinal gluons and ghost 
fields, one has to consider an effective propagator similar to those we have discussed before. 

We will first report on the result of ref . |£J , where, in the same spirit as discussed for the 
electromagnetic case, a high-temperature self-consistent QCD calculation was carried out. 

Various self energy diagrams have been taken into account in this work, including - 
compare this to the QED gluon Fock diagram for the quarks as well as polarization 

functions for the boson fields. The photon Fock diagram may be neglected in this temper- 
ature regime, since its contribution is much smaller than the gluon Fock contribution (see 
fig. |6|). Fig. [| contains a list of the diagrams, which were self-consistently summed in the 
infrared limit. Within this calculation, the quark width is obtained as 

7 9 = 0.271 gT (for N f = 2, N c = 3) . (26) 

To the same order of accuracy, and in view of the explanation in section II. B, we use the 
running coupling constant 
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FIG. 6. Contributions to the width j q of a quark. Thick dashed line: total quark width, 
calculated with the NJL model at low temperature and inspired by perturbative QCD at high 
temperature. 



9 2 {Q 2 



47T • 



12vr 



(33 - 2Nf) log Kt/Aqqj) 



(27) 



with Nf = 2 for up and down flavors, only and Aqqtj=0.2 GeV. This gives the proper match 
to the two-flavor NJL calculation towards low temperatures. Due to their larger current 
quark mass, the strange quarks will not give a large contribution at the temperatures we 
consider. We relate the mean Q 2 to the temperature by < Q 2 >~ (3T) 2 . The resulting 
quark width reads 



2.858 T 



1.65 T for T » 0.2 GeV 



log 



15.0 • T 
GeV 



(28) 



However, the question remains whether the result 7 oc gT does not contradict the calcula- 
tions in the NJL model, which yielded 7 oc aT oc e 2 T. On one hand, we may refer to eq. 
(33) of 



where 



Sg 2 T 2 
64£„ 



(29) 



Substituting for E q the minimal distance between the origin in the complex plane and the 
pole of the spectral function, i.e., E q ~ j q , we find that the self consistent •y q ~ 0.216 gT. 
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In this framework one would therefore find, that it is the replacement of m q by 7 g in the 
infrared screening, which provides a 7 9 oc gT instead of g 2 T. 

On the other hand, there are other calculations within the hard thermal loop resumma- 
tion method which give a fermion damping rate of order g 2 T [p0 |, although one may criticize 



the propagators used there because they violate the locality axiom of quantum field theory 
There, the quark width is obtained as the solution of an equation of the type []2"0| , |3"(J| , |3"7| 



gT 

g 2 m fdk fk\ g 3 T 2 ^ / g 2 T 2 l\ 







where Q(z, s, a) denotes the Lerch transcendent, a generalization of the polylogarithm func- 
tion and not expressable as a finite series of elementary functions. Numerical evaluation 
gives a value of 

7 9 w 0.82 T for T w 0.2 GeV . (31) 

For small g ~ 0.5, one finds 




7 ~ log f - I ! I - log 




(32) 



where e denotes Eulers constant. This is a non-analytical function around the point g = 0, 
but it is not of the order g 2 \og(l/g) as claimed in refs. [pO|, f2T| , j37jl (this logarithmic piece is 
also denounced in ref. ||30||). Moreover, there are hints that the logarithmic correction is due 



to a neglection of vacuum parts in the evaluated diagrams | 29| . 

The perturbative spectral width of the quark obtained in the region of the chiral phase 
transition temperature is also compatible with the high temperature NJL result, which above 
the Mott transition temperature gives a spectral width of 7 ~ T. 

Hence, starting from three completely different methods: a. the NJL model, b. effective 
QCD with generalized free fields and c. hard thermal loop resummation, we arrive at similar 
results 7 9 ~ (0.82 . . . 1.65)T in the region of the chiral phase transition. Each of these curves 
matches with the NJL result at temperatures of m 220 MeV, as can be inferred from fig. |6|. 
For reasons to be made clear later we therefore retain the parametrization of eq. Q28fl, but 
emphasize again that we do not decide at this point whether the quark damping rate is of 
order gT or g 2 T. Rather, we consider this only as a numerical parametrization which is 
supported by all available methods. 

The resulting spectral width of the quark is plotted in fig. [| together with the low tem- 
perature result obtained in the previous subsection. This temperature dependent spectral 
width of the quark therefore has two regions which we consider to be safely established 
beyond any questions of the detailed model and technique: 

{ 0.004 ■ T for T <T C (electromagnetic) 
1.0 • T for T ~ T M (strong) ^ ' 

The pronounced rise of the width around the phase transition temperature, which is carried 
over to a similar rise in the photon rates, is therefore a model independent result of finite 
temperature QCD. We identify it with the dissociation of the mesons, dominantly 7r <-» qq, 
which is connected with the phase transition temperature. 
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FIG. 7. Photon production rate R as function of the photon energy for different temperatures. 
Thick lines: Temperatures 50 - 250 MeV (dashed), our calculation using the NJL model. 
Thin line: Cutoff point Ey = 2 , y q . 



C. Results for photon production rates 

First we discuss the photon emission rates obtained from eq. ([Top below the chiral phase 
transition temperature. They are plotted in fig. ^, and we find a great similarity of the 
rate with the results plotted in fig. |2] for the purely electromagnetic case. However, at 
temperatures T > Tm ~ 212 MeV, the quark width parameter is dominated by the mesonic 
contributions, which leads to a much higher saturation scale 7 g 3> 7| m . In fig. this rise is 
also documented by the turning of the curve Ey = 2 r y q , see the thin continuous line. 

We now turn to the region of the chiral phase transition. In fig. ||, we compare the NJL 
calculations at temperature T=200 MeV (250 MeV) with the corresponding perturbative 
QCD calculations. The resulting photon production rates are different at T=200 MeV, i.e., 
very close to the chiral phase transition temperature. The reason for this is clearly, that the 
perturbative QCD calculation is no longer applicable here. 

Very good agreement is reached at a temperature of T=250 MeV, which is of course due 
to the matching values of j g at this temperature, see fig. || 

In fig. ^|, we have plotted the perturbative QCD results for temperatures T=250 and 400 
MeV, i.e., in a region where one would not trust the NJL model. The figure also contains 
a plot of the photon production rate obtained with the method of hard thermal loops, eq. 
(§). We find, that our QCD inspired calculation at the lower temperature agrees very well 
with the NJL result, while at the higher temperature and for hard photons is reproduces 
the result of the hard thermal loop calculation. For reason of this agreement we have kept 
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FIG. 8. Photon production rate R from a QGP as function of the photon energy for different 
temperatures. 

Thick lines: T=200 and 250 MeV (dashed), our calculation using eq. (28). 
Thin lines: T=200 and 250 MeV (dashed), our calculation using the NJL model. 



our parametrization for the quark spectral width above T c , eq. (|28|). 

At lower photon energies however, where according to our result the photon radiation 
is cut off due to the finite mean free path of a particle, the result of the hard thermal 
loop calculation is not usable. We therefore consider our result an extension of commonly 
accepted results to a wider range of energies. 

The leveling off of the rates at low photon energy is due to the inclusion of coherence 
(LPM) effects. For lepton pairs, this effect can be compared to the results of ref . , where 
this effect was introduced "by hand" , whereas in our calculation it is automatically included 
in the formalism. 



In fig. [10], we show the photon emission rate at three different photon energies as a 
function of temperature. Comparing the electromagnetic case (thin lines) to the model 
including the quark-meson interaction, we find a surprising result: In the region of the chiral 
phase transition, the low-energy photon production rate drops with increasing temperature. 
Eventually the radiation rate is degenerate for all energies < 2^ q (see the flat behavior of 
the curves in fig. [?]). In view of eq. ([l?]), this is understood as a dominance of the saturation 
effect over the increase of temperature. 
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Ey [GeV] 

FIG. 9. Photon production rate R from a QGP as function of the photon energy for two 
different temperatures. 

Thick solid lines: T=250 and 400 MeV, our calculation using perturbative QCD, eq. (p8| ) 
Dashed line: T=250 MeV, our calculation using the NJL model 
Thin lines: Calculation with the method of hard thermal loops, ]7| 



D. Relevance to experimental data 

We now briefly give an overview of existing or future experiments, in the order of in- 
creasing energy, which observe photons in reactions of hadronic character, and discuss the 
relevance of our rate calculations to them. 

Let us first consider very soft (for an experimentalist's scale) photons, where E 1 ~ 1 - 100 
MeV. Photons in this energy region have been measured in several experiments |2|d4,39,40 



While some of these experiments find a complete agreement of the measured photon spectra 
with the emission from hadronic sources and QED bremsstrahlung, an enhancement in the 



low p± region was observed in reactions such as K + p and n~p |40 

At the moment the discussion of experimental results is not yet conclusive, see ref. |4T 
for an excellent review of the data. However, there is one proposed explanation for such a 
soft photon excess in case it is present: It might be due to "thermal" radiation of cold drops 
of quark-gluon plasma, which would hadronize only slowly and thus have a long time to 
radiate |2|,|T|. 



While we think it premature to draw definite conclusions from the puzzling experimental 
situation, we will comment on the proposed theoretical explanation of an enhanced soft 
radiation from a cold plasma droplet. Since the emitted photons are soft, the (LPM-) effect 
of interference between successive emitters is necessarily very strong. Using our results for 
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T [GeV] 

FIG. 10. Photon production rate R from as function of temperature T. 
Photon energy 10 MeV (solid), 100 MeV (dashed) and 200 MeV (dash-dotted). 
Thin lines: purely electromagnetic width 7g m 
Thick lines, curves from left: Our calculation using the NJL model 
Thick lines, curves from right: Our calculation using perturbative QCD, eq. (pS|) 



the rate as shown in fig. |10], we therefore estimate the expected photon yield from such a 
cold plasma drop of some size R at a temperature T. From uncertainty, R > 1/p where p is 
the mean momentum of a plasma particle, and thus cold drops need to be of considerable 
size. For a drop of volume V and lifetime r at a temperature T, the differential cross section 
of photon production at a transverse momentum p± can be expressed as 

1 ' ;2Tip ± a AB R(E^T) (34) 



dp ± (0.2/ra) 

where R(E^, T) is the invariant rate and oab is the total cross section of the reaction A + B 
under study. 

As an example, let us consider the reaction K + p — > 7V where the pj_ spectrum of 
photons was measured [ 39| . At a temperature of 50 MeV, the invariant rate of photons 



of say E 1 = 10 MeV is, reading from fig. jT^, about 10 mb/GeV. Taking a drop lifetime 



r ~ R ~ 10 fm and a cross section of (Jk+ p = 16 mb results in an estimate of the photon 
yield of da/dp± ~ 3 ■ 10~ 2 mb/GeV. This is to be compared with the measured value of 
da/dp±\ eX p ~ 200 mb/GeV, which is much larger than the rate estimated from a cold plasma 
drop emitting photons. 

As we noted before, previous calculations in the cold plasma droplet picture had arrived 
at numbers which correspond to the excess over the bremsstrahlung found in some experi- 
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ments. However, these calculations do not account for the coherence (LPM-) effect which 
dramatically reduces the rate, and actually the normalization of the spectra obtained in 
these calculations was taken from experiment. We thus have to conclude, that the mecha- 
nism of "thermal" radiation from cold plasma droplets does not account for an excess of soft 
photons over the bremsstrahlung - independent of the question, whether such an excess is 
found experimentally or not. 

Now consider higher energies of the photon in the range E 1 > 100 MeV. In recent 
heavy ion experiments at ultrarelativistic energies, it is hoped to find some hints of a phase 
transition the system might, possibly partially, go through. Apart from measuring photons, 
experiments also observe lepton pairs which suffer less from coming together with a large 
background. For both electromagnetic probes, an enhancement might hint at the new phase. 



For photons, the invariant rate as shown in figs. |9| and |10| gives our result for the QGP, 
and needs to be folded with the space-time evolution of the system such as calculated in || . 
For this purpose, the invariant rate may be written in terms of the photon rapidity y and 
transverse momentum p± as 

^ da , . 

d 3 p 2irp t dp t dy 

Photons with a low virtuality can be converted into dileptons by use of the soft photon 
approximation, 

dNu a dN 

Rl+l ~ = E+E -dpidV- * "d^k ' (36) 



and improved versions thereof [ 43fl . This allows the use of the results presented in this work 
to the calculation of dilepton rates as well. Again, a space-time integration needs to be 
performed to compute the yield for a heavy ion reaction at some impact parameter, which 
is related to the measured multiplicity or total transverse energy. 

This integration over the space-time evolution of the collision was performed by use of 
the photon rates calculated from hard thermal loops |§. The system proceeds from an 
initial QGP through a mixed phase to a purely hadronic phase in the final state. Since 
the initial plasma phase is short-lived and of a similar temperature than later phases under 
the conditions studied, and the QGP does not shine very much brighter than a mixed or 
hadronic phase of the same temperature, the yield of photons from the QGP is much smaller 
than that of the other phases, typically more than an order of magnitude. As can be seen 
from fig. ^ our results in the range of photon energies E 1 > 100 MeV are similar to those 
from hard thermal loops which had been used in the analysis |§. Therefore, the conclusion 
remains, that the scarce photons from a plasma phase are overwhelmed by those from the 
later stages of the reaction, and the same is the case for dileptons. This applies to the 
application of our result to current experiments such as Ceres, WA80/98 and Helios. In 
particular, an enhancement seen in these experiments can not be accounted for by a direct 
contribution of the plasma phase, but must be of different origin (which, of course, may still 
be related to a phase transition). 



A substantial transverse flow in these collisions, for which there is some evidence |TJ| , has 
the tendency to increase the apparent photon energy and thus to increase the photon rates 
from a QGP |[45||. However, we expect that transverse flow has a similar effect on photons 
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produced from the hot hadron gas in the final stage. Therefore, this effect is probably of no 
help for the distinction of these photon sources in present measurements. 

This situation might change in favor of the quark-gluon plasma when going to higher 
energies. Here the rate of electromagnetic probes originating in the plasma rises strongly 
with increasing temperature of the QGP (see fig.|^), while the temperature at which the 
hadronic reactions occur does not change. Experiments which are under preparation at 
RHIC and LHC are planning to observe photons and dileptons and thus it is hoped that 
these experiments might see the QGP in sufficiently bright light in order to uniquely identify 
this phase. 

IV. CONCLUSIONS 

In this work, we calculated gauge invariant rates for the production of low and high 
energy photons from a hot plasma in the framework of thermal field theory with generalized 
free fields. For illustration, we studied a QED plasma and discussed possible astrophysical 
consequences. The main purpose of this work, the calculation of photon rates from a quark- 
gluon plasma, was achieved over the entire range of temperatures in a composition of two 
scenarios. For high temperature, perturbative QCD has been used, while around the chiral 
phase transition region, the nonperturbative NJL model was employed. Both results were 
found to match smoothly. This is a very satisfactory result and might be of more general 
relevance than in this particular case. For high photon energies, our result are similar to 
those obtained previously in the hard thermal loops technique, which is applicable only in 
the high temperature regime. 

It was one of the main motivations of the present work to demonstrate how meaningful 
production rates may be obtained at finite temperature for soft photons, where the coherence 
(LPM) effect plays an important role. We emphasize that the qualitative properties of the 
soft photon rates, such as the saturation effect towards low temperatures, follow from general 
physical considerations as we discussed, and are in particular independent of the particular 
model we used. In particular, the decrease of the production rate of soft photons in the 
temperature region of the phase transition is a very intriguing result, which also might have 
observable consequences. 

We pointed out that at presently reachable energies, photon production in an ultra- 
relativistic heavy ion collision is dominated by later phases of the reaction rather than an 
initially present QGP. However, at the energies of currently planned experiments, the plasma 
temperatures might be high enough to allow a direct identification of this phase. Here, the 
precise rate for the production of photons of a given energy as we calculated it is a very 
important tool for the for the interpretation of the experimental results. 

From our results one may furthermore conclude, that quantum field theory in terms of 
generalized free fields with reasonable parametrizations of spectral functions is a valuable 
method for the analysis of relativistic heavy ion collisions. The strong gap between quantum 
field theory as a formalism and its predictive power for many-body experiments, which has 
persisted for some time, is hopefully bridged by the application of thermal field theory. 
We are currently undertaking an effort to derive simple parametrizations of the photon 
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production rate that might be used as input for simulation codes. 
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